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We study black hole evaporation of near-extremal black holes in spherically reduced models 
with asymptotically Minkowskian spacetime, with the effects of the back-reaction on the geometry 
included semi-classically. The stress-energy tensor is calculated for null in-falling observers. It is 
shown that the evaporation proceeds smoothly and there are no instabilities of the outer or inner 
apparent horizon before the endpoint of evaporation. 
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The evolution of near-extremal black holes has important consequences for the resolution of the information paradox. 
Extremal black holes are potential candidates for endpoint remnants jl]|^] , and can potentially store information behind 
the double horizon || . In addition, study of their evaporation avoids the naked singularities observed in models of 
[T i • linear dilaton black hole evaporation Qj^], which are likely typical of the evolution of uncharged black holes.. 

However, it has been suggested by Jacobson, using adiabatic arguments, that the semi-classical evolution of near- 
extremal black holes may break down while still far from extremality Q . The claim is that in- falling photons created 
at the outer apparent horizon during Hawking evaporation will eventually fall through the inner horizon and either 
■ encounter a large buildup of energy behind the inner horizon, or alternatively, pile up behind the outer horizon, 
' causing it to become unstable. Thus, in either scenario, the semiclassical approximation may break down long before 
CD one would expect based on thermodynamic/statistical mechanics arguments 0. 

This claim has been explored in an analytical model § using a near-horizon Reissner-Norsdrom model originally 
proposed in [p|To|], and further elaborated in [ pd|JT^ ]. No instabilities were observed in that case. However, the 
boundary conditions there were computationally motivated rather than physically motivated. In particular there was 
| no obvious extension to asymptotically flat spacetime. 

In the present work we numerically explore semi-classical evaporation in two spherically-reduced two-dimensional 
models of asymptotically Minkowskian near-extremal black holes. We find results consistent with indicating that 
semi-classically there are no horizon instabilities during the evaporation of charged black holes. 



II. THE RN MODEL 

For our first model |L3]-[l5| , we begin with the Einstein-Hilbert action 

S = ^- J d 2 xy/=y~e- 2lt > [R + 2(V0) 2 + 2A 2 e 2 ^ . (1) 
To this, we add an electromagnetic contribution: 

Sem = -77- / d 2 x^Q 2 e 2 ^ , (2) 



2tt _ 

where Q is the charge of the black hole as it appears in the spherically symmetric Reissner-Nordstrom metric, 



The effective two dimensional action that describes (||) is given by 
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S = J d 2 x y /^e~ 2 ' t ' [R + 2(V0) 2 + 2A 2 e 20 - Q 2 e 40 ] . (4) 



In order to study the dynamics, we add a large number N of matter fields, /, 

N 



Sm — — 



^E/^^CV/,) 2 . (5) 



This form of coupling is mostly justified by the ease of the resulting solutions. However, it also corresponds to the 
bosonized version of the Callan-Rubakov modes of fermions coupled to a magnetically charged black in four dimensions 
& 

Additionally, we include the Liouville-Polyakov term 



Nh 



- / dPxy/=gRn- 1 R , (6) 



which effectively accounts for the back-reaction of the matter stress-energy tensor on the space-time geometry [|19|,[20[ . 
The full lagrangian in the conformal gauge becomes 

C = e- 2<t > (4d+d^p - 4d+<pd-4>) + X 2 e 2p - Q 2 e 2p+2cf ' - 2 K d + pd- P . (7) 

^From this point on, we choose A = 1. 



A. Equations of Motion 

Variation of ([?]) with respect to <f> gives 



(l 2 

d+d-p - d+d-4> + d+(j)d-4> + \e 2p +^ = . (8) 



A p variation gives 

" e^d+d-p- <9 + <9_0 + 2<9 + 0<9_0+ ^-(e 2 ^ - Q 2 e 4 ^) . (9) 



K ~ ■ P 2p 
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The resultant equations of motion can be written as 

d + 4>d^+\e 2p {e 2 ^~2Q 2 e^) 
d + d -P = flg^ . ( 10 ) 

(I 2 

d+d-(f> = d+d-p + d + <f>d-<f> + ^- e 2p + H . (11) 

We may note that this model possesses a singularity for e~ 2 ^ — k/2. In our solutions, this singularity will lie 
safely behind the inner horizon. It is useful to rewrite these equations with <fi — e~^. Thus d+4> = — d+(j> (f> and 
d-d + 4> = —d-d + 4> 4> + d-([)d + (j> <j>. We can rewrite (0), 



Y2 



d +( j>d_(j)+^(e 2 * -2Q 2 e^) 



d + d^ - d + <j>d^ = \e 2p ^ + 4 l 24> ^ ; , (12) 



1-fe 2 



as 



d+d-<t> Q 2 e 2p ( d+4>d-4> e 2p Q 2 e 2p 



(13) 



Hence 



4 04 \ (f> 2 4<j> 2 2</> 4 / \ 2(f> 2 



4 03 (j) 2 -K/2 
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and. 



More concisely, 



O 2 

d+d-p = d+d-<i)-d + (i)d„(j)-^-e 2p+4 ' ) ' (15) 

d+d-4> Q 2 e 2p .... 

= 7 r?r- ■ ( 16 ) 



= ^ — ' (17) 

d + d-i = -(d + d-p)i-Q^- . (is) 



B. The Static Solution 

We first solve the equations of motion in the limit of the static extremal black hole |l3| . (See also for a recent 
discussion of these static solutions, and [fl8| for analysis of the full four-dimensional semi-classical equations of motion, 
in the near-horizon limit.) We will later need these as boundary conditions for the dynamic solutions. The model 
possesses a linear dilaton vacuum solution corresponding to 

<f>=-log~(x + -x-) , (19) 
p = . (20) 



This suggests to us a reasonable choice for the radial coordinate: 



We note that 



The equations of motion, (fL7Hlq), become 



da da 1 



dx + dx~ 2 



with the constraint that 



(21) 



(22) 



(0') 2 -e 2 '(l-H!) 
P = 7 , \ * ' , (23) 



4> 2 - k/2 
~ 2 e 2p 

f/'t+OJ^. (24) 



20 U" - 2p'4>") + k (p" - {p'f + 1) = . (25) 
Linearizing about (|l9|), as a — > oo we see that 

cr 2 bp" = 1 + 2(50' - (1 + 25p) (26) 
= 250' + 26p , (27) 

and 

50" = -aSp" . (28) 

We can see that 
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6<j>" = 28 p' = -aSp" (29) 

has the solution 

M 

S P '=^, (30) 

Sfi = (31) 
a 

As a — > oo the linearized vacuum solutions are therefore 

M , s 

P = , (32) 

a 

4> = a - 2Aflogcr . (33) 



C. The Boundary Conditions 

We will study black holes perturbed by shock waves of the matter fields Our boundary conditions are chosen 
such that to the past of x + — 0, we have the static extremal solution of ((23), ( pif ) with the initial conditions at a 
chosen large value of a = determined by ([32]), ( |33| ). In practice, we make choices for Q and k, and adjust M 
so that a double horizon appears in the interior region of the solutions. Along a large value of x~ = x^, and to the 
future of x + = 0, we then impose the asymptotic linearized solutions of a black hole with mass M + AM 

*>->£^. 

5+0 = ^ ~ 2— — . (35) 

2 a; — x 



III. THE DW MODEL 

The second model we consider, the so-called DW model, derives from a general class of two-dimensional renormal- 
izable generally covariant field theories p3| , pT[ | 

C cl - ^—g{D{4>)R + G(^)(V0) 2 + H{4>)) . (36) 



We require that the potentials in (|36|) behave asymptotically like those of the linear dilaton gravity model considered 
above. That is, 

m ^m_m_ e _ 2t (37) 

as <f> — * — oo. After requiring, without loss of generality, that G(<j>) = —2D'{<f) and performing a Brans-Dicke 
transformation on the metric g = e~ 2< ^g, the Lagrangian can be rewritten as 



Cd = V-9[D{tl>)R + W(<t>)J , (38) 

where W(4>) = e 2( ^H((j)). This is the form of the lagrangian we refer to as the DW model. For these solutions, it is 
convenient to work in the conformal gauge, 

ds 2 = -e 2p dx + dx- . (39) 
Again, we add N matter fields, (||), and a back-reaction term, (||). Now 

N 

C = C c i+Y. d +fr d -fr - ^d+pd-p (40) 

i=l 

W N 
= ADd+d-p + AD'd+(j)d-(j)+ — e 2p - 2<l> - 2nd + pd- P + d+hd-h ■ (41) 

<=i 

where we have let k = N%/12. 
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A. Equations of Motion 



Variation of ( ftO| ) with respect to p and (f> yields 

4D"d+cj)d^(f> + 4D'd+d-(f> + We 2 ^ + And+d-p = , (42) 

and 

AD'd+d-p + 4D"d+<j>d-ct) - 8D'd+d-</> - 8D"d+<j>d-<j> - (W - W /2)e 2 "- 20 = (43) 

respectively. Reexpressing these in a more useful form, we obtain the equations of motion: 

D»d + <t>d-<i> + \fp-^W + Kd+d-p ( . A . 
d+d-cj) = ^ , (44) 

^ = ^t^ ■ (45) 

As discussed in pi] , we must choose D such that D'{4>) < for all <^> in order to avoid a singularity of the Brans-Dickc 
transformation. Also, 

g ~ J d<pW{4>)D'{<p) (46) 

has zeroes along a, the radial coordinate, equal to the number of zeroes of W((p)D'(^>) plus one. Therefore, for a 
charged black hole, W ((f>) D' ((f>) must possess a single zero, for which <j){cr) will correspond to the extremal radius. In 
our solutions, we investigate the choice where D = e~ 2 ^ — 7 2 e 2< ^ and W — 4 — p 2 e A ^ . p 2 , as the coefficient of e 2 ^, 
behaves as the charge of the black hole. This becomes clear through comparison with (||), from the Reissner-Nordstrom 
model. 

B. The Boundary Conditions 

The boundary conditions are set so that to the past of Xq = 1 the solution corresponds to the static, extremal 
solution of the equations of motion, while along x~ = — oo, and to the future of Xq = 1, the boundary conditions are 

(j> = p = -^log (M + AM(4 - x + ) - x + x~) , (47) 

_ _ x" + AM 

+P ~ 2(M + AMx£ - x+x- - AMx+) ' ^ ' 

This is the same boundary condition imposed in In practice, the static solutions are solved with the static form of 
([l||) and (|4^), imposed at some large value of x~ = x~ Then, M is fine-tuned so that the there is a double horizon 
static solution to the equations of motion, (|44|) and (|45|). This solution in turn becomes the boundary condition at 
Xq = 1 for the full dynamical solutions. 



IV. AFFINE COORDINATES 



We would like to answer questions regarding black hole stability during the evaporation process. An appropriate 

indicator would be the behavior of the stress-energy as observed along a null affincly parameterized geodesic. The 
Christoffel symbols for the metric ([39|) are 

T+ + = 2d +P , (49) 

r:_ = 2d-p . (50) 

^From the geodesic equation for x + , we obtain 

d 2 x + _„ dx + dx + „ 

dx^ +2d+P dxTdxT = °> < 51 > 
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with a similar equation for x . x + is an affinely parameterized null geodesic. If x is held fixed at a; = x , 



Then (Ell) reduces to 



Thus, 



which is then solved to give 



d 2 x + dp dx+ _ 

dx+t + 2 dx+dx+~ U ■ [M) 



^ (54) 



^=C(x^)e 2 ^ + ^) , (55) 

where C(x~) is a constant of integration which we will soon determine. Similarly, at fixed x + = Xq , 

dx~ 
dx~ 

Thus we make calculations of T, as defined by 



C(x+)e 2p ( x o' x > . (56) 



T±±=[^) T ±± --^T ±± . (57) 

This gives us an observable measure of the energy flux that is independent of the choice in coordinates. 

We choose the normalization constants, C(xjT), C(xq) in ( |55| ) and (|56|), so the affine coordinates coincide with 
asymptotically Minkowskian coordinates far from the black hole. This normalization will not be necessary in the RN 
case, because there the coordinates are already asymptotically Minkowskian. For the DW model, we choose C(Xq), 
C(xq) with the requirement that along some large chosen value of x + , termed x + — x+ , and along a similarly chosen 

ds 2 = -dx+dx- , (58) 

so that, 

e2p dx+ dx_ = 1 
dx + dx~ 

Given that asymptotically e~ 2p = M — x + x~ , then this amounts to requiring that at x~ = Xq 



dx+ »-2o(x+.x7)VM- 



= e-W* + <** > v 111 7_ 3 ° , (60) 
u,it ■ M — xXqXq 

and at x + = Xq 

^1(4) = e -2 P (4, a -)V ^-^^° . (6i) 

Using these we obtain meaningfully normalized null affine coordinates with which to calculate the stress-energy tensors 

of m). 
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V. RESULTS 



Careful analysis limits the range of values of our variables, (7, fj, and K, in the DW model, Q and k in the RN) for 
which numerical errors are sufficiently small to permit unambiguous statements. We compared our results against the 
classical solutions to gauge the buildup of numerical error. These numerical constraints prevent us from making the 
back-reaction arbitrarily weak. A useful measure of the strength of the back-reaction is the ratio of the time-scale of 
the quantum evaporation to the light-crossing time of the black hole. These time-scales may be extracted by inspection 
of the equations of motion, with the results that tu g ht/t e vap — k/Q 2 for the RN model, and tu g ht/t e vap = for 
the DW model. We quote the results for the regime when this ratio is of order one, where numerical errors are 
negilible. The qualitative behavior of the solutions does not change as n is made smaller, as far as we have been able 
to determine. Of course, the semi-classical approximation to the evaporation is no less valid when we do not impose 
the weak back-reaction condition. 
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FIG. 1. DW model: The outer horizon is shown as it recedes to meet the inner horizon, in turn moving out, at <fio, the 
extremal radius. x + and x~ are the x and y axes respectively. [7 = 8; fj, = 15; k = 10; AM = 1.5] 

A typical contour plot is shown for each model in FIG. [l] and FIG. ||. As a shock mass is sent in, the two apparent 
horizons, given by the zeroes of d+(j>, split apart, with the outer apparent horizon moving past (f>o, while the inner 
horizon moves behind this radius. Immediately afterwards, the outer horizon begins to recede and meets the inner 
horizon back at the extremal radius <po (the value of <fr at the position of the horizon of the initial extremal black 
hole). 

It is useful to regard the separation of the horizons along the x~ direction as a measure of the excitation energy 
of the extremal black hole (this statement can be made more precise in the adiabatic approximation Jl5| , p^| ). This 
indicates the semi-classical solutions break down before the meeting point of the apparent horizons since here the 
energy of an emitted quanta inevitably becomes comparable to the energy of excitation above the extremal ground 
state . Evolving for points in the causal future of this endpoint would no longer be consistent with the semi-classical 
approximation. 

As discussed in the Ricci scalar curvature calculated along a given contour of <f> (see FIG. || and FIG. ||) 
demonstrates a return to the extremal quantity as evaporation proceeds, strongly indicating a smooth evaporation 
process. This is to say that an observer at a fixed radius outside of the black hole sees the same physical environment 
towards the endpoint of evaporation as was seen before the shock mass was introduced. 

For the DW model, we also evaluate T±± (as shown in FIG. || and FIG. ||), calculated from ( |57| ) along the fixed 
4> contour in FIG. [l|. The behavior of these quantities is similarly consistent with a smooth evaporation process, 
returning the near-extremal black hole back to the extremal ground state. These fluxes are non-zero initially, but 
approach zero during evaporation. In the RN case, however, the flux perceived by an observer at a fixed radius outside 
of the black hole does not go to zero. This is evident from inspecting the behavior of T__, shown in FIG. ||, whose 
absolute value does not decrease monotonically, indeed eventually rising again. This is a not an entirely unexpected 
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FIG. 2. RN model: The outer horizon is shown as it recedes to meet the inner horizon, in turn moving out, at <fio, the 
extremal radius. x + and x~ are the x and y axes respectively. The additional contour lying above these is a remnant of (j>o. 
No particular meaning has been attributed to this as it lies in a causally inaccessible region. [Q = \/60; k = 40; AM = .2] 
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FIG. 3. DW model: The Ricci scalar curvature, R, calculated along the contour <fi = —1.3 throughout the evaporation 
process. The asymptotic value to which R settles coincides with its extremal value. [7 = 8; // = 15; k = 10; AM — 1.5] 
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FIG. 4. RN model: The Ricci scalar curvature, R, calculated along the contour <f> — 15 throughout the evaporation process. 
The asymptotic value to which R settles coincides with its extremal value. [Q = k = 40; AM — .2] 
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FIG. 5. DW model: T ++ calculated along the contour <f> — —1.3 varies smoothly throughout the evaporation pro- 
cess. The data suggests that the flux T ++ approaches zero, or a small constant, as the shock mass is evaporated away. 
[7 = 8; n = 15; k = 10; AM = 1.5] 
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FIG. 6. DW model: T__ calculated along the contour <j> — —1.3 varies smoothly throughout the evaporation pro- 
cess. As the endpoint of evaporation is reached, the data suggests that this quantity approaches some constant value. 
[7 = 8; fj. = 15; k = 10; AM = 1.5] 




FIG. 7. RN model: T ++ calculated along the contour <f> — 15 throughout the evaporation process. After the shock mass is 
introduced at Xq , the flux T ++ , which is negative, effectively approaches zero. [Q — \/60; k — 40; M — 5.641; AM — .2] 
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FIG. 8. RN model: T — calculated along the contour cj> = 15 throughout the evaporation process. 
[Q = V60; k = 40; AM = .2] 
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FIG. 9. DW model: The Ricci scalar curvature, R, evaluated at points along the outer horizon, the inner horizon, and 
4>o, overlapped. The absolute value of the curvature at the outer horizon is increasing, while decreasing at the inner horizon. 
When the horizons meet again at the extremal radius, <f>o, the curvature is at its extremal value, indicating that the horizon 
environment has returned to its original state. [7 = 8; fj, = 15; n = 10; AM = 1.5] 
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FIG. 10. RN model: The Ricci scalar curvature, R, evaluated at points along the outer horizon, the inner horizon, and 
4>o, overlapped. The absolute value of the curvature at the outer horizon is increasing, while decreasing at the inner horizon. 
When the horizons meet again at the extremal radius, <f>o, the curvature is at its extremal value, indicating that the horizon 
environment has returned to its original state. [Q = \/60; k = 40; AM = .2] 
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result, since the affine coordinates are not true asymptotically Minkowskian coordinates, so the flux need not vanish 
at infinity. T++ (FIG. @), however, does approach a zero value with increasing x . 

In order to directly test for instabilities in the inner or outer apparent horizon we also inspect the behavior of 
the stress-energy momentum and the Ricci scalar curvature for the contours of d+4> — along the outer and inner 
apparent horizons in order to see if an instability arises. FIG. [)] and FIG. [lO] appear to indicate that the scalar 
curvature evolves smoothly throughout evaporation and returns to its extremal value. Meanwhile, FIGs. [Tl] - [l8| 
demonstrate that for both models the absolute flux through the horizons decreases during the evaporation, consistent 
with a decreasing black hole mass. More importantly, there is no indication of a buildup of energy behind either the 
inner or outer apparent horizon. The smooth variation of T ++ and T__ in these regions suggests a return, through 
evaporation, to the extremal ground state after the initial excitation. 
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FIG. 11. DW model: T++ evaluated at points along the outer horizon, shown as it varies smoothly for increasing x + . 
[7 = 8; n = 15; k = 10; AM = 1.5] 
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FIG. 12. DW model: T__ evaluated at points along the outer horizon, shown as it varies smoothly for increasing x + . 
[7 = 8; n = 15; k = 10; AM = 1.5] 



VI. DISCUSSION 



The pictures emerging from both the DW and RN case were consistent. After the injection of shock matter, the 
near-extremal black hole evaporates and returns to the extremal ground state. The outer and inner horizon experience 
no distressing energy instabilities and come back together, rejoining at the extremal radius. Furthermore, studies of 
the Ricci scalar and the stress-energy tensors calculated along null affine geodesies, indicate that for observers outside 
the black hole, the environment returns to its original state. 

The original motivation for studying the evaporation of charged black holes was to shed light on the black hole 
information puzzle. Quantum mechanics necessarily requires unitarity. However, semi-classically, the emitted radi- 
ation from black holes possesses a thermal distribution. This leads to the information puzzle: what happens to the 
information contained in the ingoing matter once the mass of the black hole is radiated away? 

While it may be that physics in the real universe is simply not unitary ]23| , |24| ] , other resolutions have nevertheless 
been suggested. One such potential resolution of the paradox is that the information of the ingoing matter becomes 
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FIG. 13. DW model: T ++ evaluated at points along the inner horizon, shown as its value decreases smoothly for increasing 
x + . [7 = 8; n = 15; k = 10; AM = 1.5] 
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FIG. 14. DW model: T — evaluated at points along the inner horizon, shown as it varies for increasing x + . This flux goes 
to zero as the shock mass evaporates away. [7 = 8; n = 15; n = 10; AM — 1.5] 
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FIG. 15. RN model: T__ evaluated at points along the outer horizon, shown as it varies smoothly for increasing x 
[Q = \/60; n = 40; AM = .2] 
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FIG. 16. RN model: T ++ evaluated at points along the outer horizon, shown as it varies smoothly for increasing x + . The 
T++ flux through the outer horizon appears to vanish as evaporation proceeds. [Q = k = 40; AM — .2] 
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FIG. 17. RN model: T++ evaluated at points along the inner horizon, shown as it varies for increasing x . This flux 
decreases in absolute value throughout the evaporation process. [Q = a/60; k = 40; AM = .2] 
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FIG. 18. RN model: T__ evaluated at points along the inner horizon, shown as it varies smoothly for increasing x + . 
[Q = n = 40; AM = .2] 
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stored in stable remnants that remain at the end of evaporation. These remnants would apparently need to contain 
an arbitrarily large amount of information despite possessing a Planck-sized mass |2(if| . In this regard, extremal black 
holes hold some promise; behind the charged black hole horizon lies an infinite throat region j|] that can apparently 
store an arbitrary amount of information. This fact is also manifested by the infinite degeneracy of the extremal 
ground state p| at the semi-classical level. Thus extremal black holes represent the testing ground for the remnant 
resolution to the information paradox §],§]. There is some debate as to whether or not the pair production of such 
charged black holes occurs at a finite rate p|,[i|p7|]2[ , given the infinite degeneracy of the ground state. However, this 
infinite degeneracy is likely a result of the semiclassical averaging over the quantum geometries. This should not occur 
in a full quantum treatment. For example, D-brane calculations demonstrate that summing over microstates in fact 
yields a finite Bekenstein-Hawking entropy |?8). Matrix theory also indicates a similar result |2fi| . This may lay to 
rest infinite production concerns. Our studies did not reveal an instability in the evolution of excited extremal black 
holes. Thus as it stands, we cannot yet rule them out as the ultimate devices for information storage. 

Another suggested resolution to the information puzzle proposes that higher order quantum corrections would reveal 
that the outgoing radiation is not strictly thermal, carrying off information about the ingoing radiation in the form of 
quantum correlations with the ingoing matter. This, however, while at first seeming to be a most agreeable resolution, 
presents other theoretical difficulties. Once matter disappears behind the horizon, it becomes causally dissociated 
from the outside region. Therefore such cross-horizon correlations would necessarily violate causality, unless some 
sort of information stripping were to occur at the horizon itself. At the semi-classical level, however, we do not expect 
anything special to occur at the horizon in the frame of reference of the in- falling observer. This issue can be resolved 
by appealing to a "stretched horizon" membrane picture |30| , based on the principle of black hole complementarity. 
This principle asserts that matter which has fallen past the event horizon and the Hawking radiation emitted from 
there are not different objects; they are equivalent and complementary descriptions from reference frames which can 
be related to each other by large Lorentz boosts. Since there is no communication to outside of the black hole once 
an observer has passed in, there is no logical contradiction by prescribing that unitarity remains valid entirely outside 
of the horizon. Furthermore, the contribution of non-local terms arising from string theory |n],[32| may account for 
cross-horizon correlations which will nevertheless be undetectable by any single observer. That is, if we do not require 
a local field theory, except at low energy, it becomes possible to construct a fully self-consistent causal theory where 
field operators behind and in front of the horizon will not commute [ j33| . String theory and the existence of quantum 
correlations between the ingoing and outgoing matter then ultimately promises to contain the key for resolving the 
information puzzle. 

The models we have considered here suggest that at the semi-classical level, upon entering the black hole, information 
is lost to the outside world. However, while it becomes inaccessible, the information enters the infinite asymptotic 
throat region behind the horizon and is not actually destroyed. While this picture permits remnants to contain the 
"lost" information, it still remains possible that the non-local string theory terms do indeed represent the real universe 
and the information does return to the outside region. Our success lies in showing that the semi-classical treatment 
of charged black hole evaporation does not breakdown sooner than would be expected from [Q, despite heuristic 
arguments of Jacobson ^] to the contrary. Results presented here are also consistent with analytical results found in 
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